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Abstract: The central exclusive production (CEP) of heavy resonance states that sub- 
sequently decay into meson pairs, MM, is an important signature for such processes at 
hadron colliders. However there is a potentially important background from the direct 
QCD production of meson pairs, as mediated for example by the exclusive gg — )• MM 
hard scattering subprocess. This is in fact an interesting process in its own right, testing 
novel aspects of perturbative QCD technology. We explicitly calculate the gg — t- MM 
helicity amplitudes for different meson states within the hard exclusive formalism, and 
comment on the application of MHV techniques to the calculation. Using these results, we 
describe how meson pair CEP can be calculated in the perturbative regime, and present 
some sample numerical predictions for a variety of final states. We also briefly consider the 
dominant non-perturbative contributions, which are expected to be important when the 
meson transverse momentum is small. 
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1. Introduction 



Central exclusive production (CEP) processes of the type 



pp{p) -^p + X + pip) , 



(1.1) 



can significantly extend the physics programme at high energy hadron colliders. Here 
X represents a system of invariant mass Mx, and the '+' signs denote the presence of 
large rapidity gaps. Such reactions provide a very promising way to investigate both 
QCD dynamics and new physics in hadron collisions, and consequently they have been 
widely discussed in the literature (see [1-5] for reviews). Of special interest is the study 
of resonance states, from 'old' SM mesons to BSM Higgs bosons (see for example [6-11] 
and references therein). An attractive advantage of CEP reactions is that they provide an 
especially clean environment in which to test the nature and the quantum numbers of the 
centrally produced state X. 

Recall (see [1,12,13]) that in exclusive processes, as shown in Fig. 1, the incoming gg 
state satisfies special selection rules in the limit of forward outgoing protons, namely it has 
Jz = 0, where Jz is the projection of the total gg angular momentum on the beam axis, and 
positive C and P parity. Hence only a subset of the helicity amplitudes for the gg ^ X 
sub-process contributes. The CEP mechanism therefore provides a unique possibility to 
test the polarization structure of the gg X reaction. 

Even without forward proton spectrometers, double diffractive processes of the type 



with large rapidity gaps separating the centrally produced system X from the products, 
Y and Z, of the proton (antiproton) dissociation, are still of considerable interest. If the 
incoming protons are dissociated into low masses {My,z ^ 2 GeV), then these reactions 
continue to exhibit many of the attractive properties of CEP, while in the case of high- 
mass dissociation the Jz = selection rule, characteristic of exclusive production, no longer 
holdsi. 

These double diffractive processes can be measured using forward rapidity gap triggers, 
with the help of simple scintillation (forward shower) counters (FSCs) [14]. Such a strategy 
was successfully implemented at the Tevatron by the CDF collaboration [5,15-18], where 
the rapidity gap trigger was used to veto on particles with pseudorapidity \r]\ < 7.5 on 
each side of the central system. CDF have published a search for 77 CEP [16] with 
^t{i) > 5 GeV. This process (together with the CEP of charmonia, the observation 
of which was reported in [17]), can serve as a 'standard candle' reaction with which we 
can check the predictions for new physics CEP at the LHC [8, 19]. Based on the results 
in [19], the observed CDF 77 CEP cross section should be small, corresponding to 0.81q5 
events; experimentally [16], three candidate events were observed. Subsequently the Et{'j) 
threshold has been decreased to ~ 2.5 GeV, and many more candidate events have been 

^See footnote 3 after (2.8). 



pp(p) -^Y + X + Z , 



(1.2) 
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observed [5,20]. In this context, a good quantitative theoretical understanding of the CEP 
background is crucial, since one or both of the photons from tt^ 77 decay can 
mimic the 'prompt' photons from gg — )• 77 CEP. 

As discussed in [8,21,22], the observation of XcO CEP via two-body decay channels is 
of special interest for both studying the dynamics of heavy quarkonia and for testing the 
QCD framework of CEP, with the Xc — ^ tttt decay being a promising example. We recall 
that these channels, especially vrvr, K^K~ and pp, arc ideally suited for spin-parity analysis 
of the Xc states. In particular, the fact that the Xc(i.2) two-body branching ratios arc in 
general of the same size or smaller (or even absent for the Xci) than for the XcO^ ensures 
that the = selection rule is fully active, see [12,21,22] for more details. However, in 
this case we may expect a sizeable background resulting from the direct QCD production 
process 

PP{P) p + hh+ p{p) (1.3) 

with h = IT, K,p. Such a non-resonant contribution should therefore be carefully evaluated. 
We recall that the existing CDF measurement of Xc CEP [17] is based on the detection of 
the decay Xc ^ J /'ijj + ^. While the overall observed rate is in agreement with theoretical 
expectations [21,22], the M(J/V'+7) resolution, due to the low photon energy (~ 200 MeV), 
does not allow a separation of the different XcJ states. Therefore, although the XcO should 
be the dominantly (CEP) produced Xc state, the higher spin Xci and Xc2 states could 
give a sizeable contribution to the observed J/tp + ^ signal, because of their much higher 
branching fractions to J/ip +'j [22,23]. 

While the Tevatron CEP data do not in general make use of forward proton detectors, 
relying instead on large rapidity gap triggers, a new area of experimental studies of CEP 
with to^^ed forward protons is now being explored by the STAR Collaboration at RHIC [24], 
which has the capability to trigger on and to measure forward protons. This provides an 
excellent means to extend the physics reach in studying CEP processes. 

Currently at the LHC, without forward proton taggers, the existing detectors (ALICE, 
ATLAS, CMS and LHCb) were not designed for studying such standard candle processes 
as charmonia (e.g. Xc-) or 77 CEP, because they lack the forward coverage necessary to 
measure large rapidity gaps. Selecting only events with a low multiplicity in the central 
detector may not be enough to study these reactions since such events include processes 
with proton dissociation, as in (1.2) However, as mentioned above, the addition of FSCs 
could provide sufficient rapidity coverage and would allow the exclusion of events with high- 
mass and a large fraction of events with low-mass diffractive dissociation. Thus events with 
a 'veto' FSC trigger (assuming that the number of interactions per bunch crossing is not too 
large) may be interpreted as arising from CEP. In fact such counters have been proposed for 
CMS [14,25] and are currently being installed. This will open up a wide physics programme 
at CMS based on vetoes using the FSC, ZDC, CASTOR and HF detectors, together with 
requiring some minimal activity in the central region [5,14,25,26]. In particular, CEP 
events of the type (1.1) can be selected without detecting forward protons. 

Finally, we note that by employing FSCs a rich CEP physics programme could also be 
realised with the LHCb experiment [27]. The excellent particle identification of the LHCb 
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detector and the high momentum resolution are especiahy beneficial for measurements of 
low multiplicity final states^. A promising study of low central mass CEP events is also 
ongoing at ALICE [29], using additional scintillator detectors placed on both sides of the 
central barrel, which allows tagging of double rapidity gap events. 

Studies of meson pair CEP would also present a new test of the perturbative formalism, 
with all its non-trivial ingredients, from the structure of the hard sub-processes to the 
incorporation of rescattering effects of the participating particles. Of particular interest is 
r]r] and rj'r]' CEP, which could allow a probe of the gluonic structure of 77, 77' mesons, the 
available theoretical information on which is still quite scarce. Though being undoubtedly 
of interest in their own right, the two-meson non-resonant exclusive processes (1.3) are of 
essential practical importance as the main backgrounds to Xc ~^ KK CEP, and, in the 
case of vrO^o, to 77 CEP. This in turn requires a very detailed theoretical knowledge of the 
properties of the process (1.3). 

TTTT CEP, mediated by Pomeron-Pomeron fusion, has been the subject of theoretical 
studies within a Regge-pole framework since the mid-1970s (see, for example [30, 31] for 
early references and [32,33] for more recent ones). There have also been a variety of exper- 
imental results on low-mass meson pair CEP, in particular from the CERN ISR, with the 
measured cross sections in broad agreement with the expectations of Regge phenomenol- 
ogy (sec [4,5] for reviews). As discussed in [33], at comparatively large meson transverse 
momenta, p±, CEP should be dominated by the perturbative 2-gluon exchange mecha- 
nism of Fig. 1, where the gg — > tttt, K'^K~ coupling can be calculated using the formalism 
of [34,35]. In the kinematic regime relevant for Xc CEP (M^^r ~ M^, p±{Tr) ~ both 
non-perturbative and perturbative mechanisms could in principle contribute to the overall 
rate, and this issue requires careful investigation. Of special interest is the transition re- 
gion between these two contributions, which is very sensitive to the behaviour of the meson 
form factor FM{t). The potential tttt CEP continuum background to XcO was also recently 
considered purely within the framework of Regge theory by [36] (see also [37]) - we will 
briefly comment on and compare our results in Section 6. 

In this paper we perform a detailed study of the CEP of meson pairs, paying special 
attention to the perturbative regime. First, in Section 2, we review the CEP calculation 
formalism for processes (1.1). We then in Section 3 consider meson pair, MM, production, 
in particular the perturbative QCD calculation of the gg — t- MM process, thus generalising 

^Recently the LHCb collaboration have reported encouraging preliminary results on possible CEP in the 
Xc(o,i,2) — > J/tp + 7 channel, where in order to justify the exclusivity vetoing was imposed on additional 
activity in the rapidity region 1.9 < < 4.9, with some sensitivity to charged particles in the backwards 
region —4 < r] < —1.5 [28]. While the Xc(n.i) production data are in good agreement (within theoretical 
and experimental uncertainties) with our predictions for CEP [8, 22] , the observed x<;2 rate is somewhat 
higher. However, as discussed above the observed data will include events with proton dissociation, as 
in (1.2). Such an inclusive contamination, where the momentum p± transferred through the t-channel gg 
system (see footnote 3) may be rather large, will tend to preferential^ increase the higher spin Xc{i,2) yield. 
Thus the xcz cross section, which is proportional to p^, could be particularly enhanced. This may therefore 
account for some or all of the observed disagreement between the LIfCb Xc2 measurement and the CEP 
theory. Although a data-driven attempt to subtract the inclusive contribution has been made, it is not clear 
that this will account for all of the background: as discussed above, the addition of FSCs on both sides of 
the LHCb experiment would allow a more efficient veto on such inelastic events [27], and would therefore 
greatly clarify the situation. 
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the existing 77 — >• MM results. Unlike the 77 case, the gg — )■ MM cross section exhibits 
an interesting (QCD) radiation zero. We also show how MHV techniques can be used 
to calculate the underlying gg — )• qqqq helicity amplitudes. In Section 4 we address the 
calculation of r]rj, r/r/' and r]'r]' production, and in Section 5 we calculate and compare the 
subprocess and hadron-level cross sections for various final states. Non-perturbative meson 
pair CEP is considered in Section 6, and a secondary perturbative mechanism is described 
and compared with our standard mechanism in Section 7. Finally, Section 8 contains our 
conclusions and outlook for further work. 



2. Central exclusive production 

The formalism used to calculate the perturbative CEP cross section is explained in detail 
elsewhere [6,8,13,38] and so we will only review the relevant aspects here. The amplitude 
is described by the diagram shown in Fig. 1, where the hard subprocess gg ^ X \s initiated 
by gluon-gluon fusion and the second t-channel gluon is needed to screen the colour flow 
across the rapidity gap intervals. We can write the 'bare' amplitude in the factorised 
form [1,19,21,39] 

T = TT^ f -2— ^ ^ fg{xl,x'„Ql^Ji''■M)fg{x2,x'^,Ql^i^;h) , (2.1) 

J Qi(Q± -Pu) (Q-L + P2x) 

where the /^'s in (2.1) are the skewed unintegrated gluon densities of the proton: in the 
kinematic region relevant to CEP, they are given in terms of the conventional (integrated) 
densities g{x,Qf). ti is the 4-momentum transfer to proton i and ji is the hard scale of the 
process, taken typically to be of the order of the mass of the produced state: as in [21,39], 
we use ^ = Mx/2 in what follows. The t-dependence of the /^'s is isolated in a proton 
form factor, which we take to have the phenomenological form F]\j{t) = exp(6t/2), with 
b = 4GeV~^. The M. is the colour- averaged, normalised sub-amplitude for the gg ^ X 




Figure 1: The perturbative mechanism for the exclusive process pp ^ p + X + p, with the eikonal 
and enhanced survival factors shown symbolically. 
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process 

M = - 

'X^'C - a,b 



Here a and b are colour indices, Mx is the central object mass, V^^! represents the gg ^ X 
vertex and qi_^ are the transverse momenta of the incoming gluons, given by 

9ix =Q±- Pi± , 92x = -Q± - P2_L , (2.3) 

where Q± is the momentum transferred round the gluon loop while pij^ are the transverse 
momenta of the outgoing protons. Only one transverse momentum scale is taken into 
account in (2.1) by the prescription 

Qi = min{Q_L, |(Q± - Pi_l)|} , 

Q2 = min{Qx,|(Qx + P2x)|}- (2-4) 
The longitudinal momentum fractions carried by the gluons satisfy 

where x' is the momentum fraction of the second t-channel gluon. The differential cross 
section at X rapidity yx is then given by 

/c-2 V /■j2„ ,2„ |T^(Plx,P2x))Pc.2 



/c2 \ f j2„ k(P1x'P2x)) c2 / „ \ /cy (i\ 

= {Senh) / a Plxd P2x S'eik(Plx>P2x) > (2-6) 



where T is given by (2.1) and S'gjj^ is the 'eikonal' survival factor, calculated using a gen- 
eralisation of the 'two-channel eikonal' model for the elastic pp amplitude (see [40] and 
references therein for details). 



Besides the effect of eikonal screening S'eik, there is some suppression caused by the 
rescatterings of the intermediate partons (inside the unintegrated gluon distribution, fg). 
This effect is described by the so-called enhanced Reggeon diagrams and usually denoted 
as Sgj^h, see Fig. 1. The value of ^^^^j^ depends mainly on the transverse momentum of 
the corresponding partons, that is on the argument of fg{x,x',Q^,iJ,^;t) in (2.1), and 
depends only weakly on the p± of the outgoing protons (which formally enters only at 
NLO). While S'pjjj was previously calculated (sec [8,22]) using the formalism of [41], we 
now use a newer version of this model [42] which includes the continuous dependence on Qf 
and not only three 'Pomeron components' with different 'mean' Qi. We therefore include 
the Senh factor inside the integral (2.1), with {S'^^-^) being its average value integrated over 

Q±- 



If we consider the exact limit of forward outgoing protons, pij_ = 0, then we find that 
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after the integration (2.2) reduces to [1] 

M oc ql^qi^V,, ^ \qI{V++ + F--) ~ J2 ^^'^'^x^x, , (2.7) 

Ai,A2 

where A^i 2) are the gluon helicities in the gg rest frame. The only contributing hehcity 
ampUtudes are therefore those for which the gg system is in a J2 = state, where the 2;-axis 
is defined by the gg direction which, up to corrections of order ~ q^/M^, is aUgned with 
the beam axis. In general, the outgoing protons can pick up a small p±, but large values 
are strongly suppressed by the proton form factor, and so the production of states with 
non- J2 = quantum numbers is correspondingly suppressed (see [8, 22] for an example of 
this in the case of X(c,b) ^^'^ ''l{c,b) CEP). In particular, we find roughly that 

\Tm = 2)\^ {plf 

|T(J, = 0)|2 ~ (Q^)2 ' 

which is of order ~ 1/50 — 1/100, depending on the central object mass and cms energy -s/s. 
We shall see that this 'J^ = selection rule' [6,12,13] will have important consequences for 
the case of meson pair CEP^. Finally, we note that in (2.7) the incoming gluon (transverse) 
helicities are averaged over at the amplitude level: this result is in complete contrast to 
a standard inclusive production process where the cross section is averaged over all gluon 
helicities. Eq. (2.7) can be readily generalised to the case of non-J^ = gluons which 
occurs away from the forward proton limit, see in particular Section 4.1 (Eq. (41)) of [8], 
which we will make use of throughout to calculate the MM CEP amplitudes from the 
corresponding gg MM helicity amplitudes. 

3. The gg — )■ MM perturbative process 

3.1 Background: large angle 77 — > MM meson pair production 

Before considering the gg — t- MM process relevant to CEP, we begin by summarising the 
case of 77 —7- MM, which has been considered previously in the literature; the formalism 

^In the perturbative domain the hierarchy p± < Q± is provided by the mass inequahty Mx ^ My,z, 
that is the mass of the centrally produced state, Mx, is large in comparison with the mass of the proton or, 
in the case of the process (1.2), the states. My, Mz, originating from the proton dissociation. Recall that 
normally in the hard process a large momentum transfer, p±, (this is the momentum transferred through the 
t-channel gg system, similar to that entering in (2.3) for the purely exclusive process (1.1)) leads to proton 
dissociation into higher mass states Y, X. Thus the probability to observe a system of low mass My,z 
when there is a large p± transfer is small. The dominant contribution in the low mass, My,z, region comes 
from nucleon resonances where large values of p± are suppressed by the nuclcon, or N ^ N* transition, 
form factors. On the other hand, the Sudakov form factor hidden in the unintegrated gluon densities, fg, 
strongly suppresses the contribution from low Q± <S Mx, leading to a rather large average Q± in the 
integral (2.1). Moreover, there are phenomcnological indications (such as the small value of the slope, a'p, 
of the Pomcron trajectory, the success of the additive quark model, a{np) ~ (2/3)a{pp), etc.) that the 
radius of the 'soft' Pomcron is small and that it should have a qualitatively similar structure to the 'hard' 
(QCD) Pomeron (see [42,43] for more details). Thus even for a relatively small Mx we may consider the 
diagram shown in Fig. 1 to have a relatively large average Q±. Thanks to the fourth power {p±/Q±)* in 
(2.8), already for Q± ~ 2p± the Jz = selection rule becomes quite effective, with the | Jz| = 2 contribution 
being suppressed by a factor of ~ 16. 



(2.8) 
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used to describe this process can then be readily applied to the gg case. The leading order 
contributions to 77 MM were first calculated in [34] (see also [35,44]), where M(M) is 
a flavour nonsinglet meson(antimeson)^. They can be written in the form 

Mxy{s,9)= f dxdyct>M{x)^{y)Txx'{x,y-s.e) . (3.1) 
^0 

where s is the MM invariant mass, x. y are the meson momentum fractions carried by the 
quarks, A, A' are the photon helicities and 6 is the scattering angle in the 77 cms frame. 
Txx' is the hard scattering amplitude for the parton level process 77 — t- qqqq, where each 
(massless) qq pair is collinear and has the appropriate colour, spin, and flavour content 
projected out to form the parent meson. In the meson rest frame, the relative motion 
of the quark and antiquark is small: thus for a meson produced with large momentum, 
\k\, we can neglect the transverse component of the quark momentum, q, with respect to 
k, and simply write q = xk in the calculation of Txy- <P{x) is the meson wavefunction, 
representing the probability amplitude of finding a valence parton in the meson carrying a 
longitudinal momentum fraction x of the meson's momentum, integrated up to the scale 
Q over the quark transverse momentum (with respect to pion momentum k). 

A representative diagram contributing to Txx' is shown in Fig. 2, where the additional 
hard gluon is required to supply the momentum transfer between the quark lines necessary 
for large angle scattering. There are four basic LO Feynman diagrams which contribute 
(the rest being given by simple permutations of the fermion lines and incoming photons), 
and after an explicit calculation we flnd^ 



Airas 327rQ! 



^ s x(l — x)y{l — y) 
+_ ^ +_ ^ „ 47ra5 327ra 
1^^ ^77 . x{\-x)y{l-y) 



(ei - 62)^0"! 



(3.2) 



(ef - el){x-y) 



1 — cos^ 9 

(ei - 62)^(0 - 1) _ eie2a(j/(l - y) + x{l - x)) 
1 — cos^ 9 a? — IP' cos^ 9 

(3.3) 



consistent with the results of [34], up to a universal '— ' sign, ei, 62 are the quark charges 
(i.e. the mesons have charges ±(ei — 62)) and 

a={l~x){l-y)+xy , (3.4) 
h={l-x){l-y)-xy . (3.5) 

A few comments are in order about the meson wavefunction ^{x). Firstly, as is usual for 
the factorisation of long- and short-distance physics, (f){x) must evolve with in order 
for the physical predictions to remain factorisation scale invariant (see for example [46] for 
a detailed discussion). In particular, ^(a;) = (p{x,Q'^) obeys an evolution equation of the 



*As pointed out in [34], see also [45], for the case of mesons with flavor-singlet Fock states there is also 

a contribution coming from the LO two-gluon component, see Section 4. 

^We have made use of the FORM symbolic manipulation programme throughout. 
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Figure 2: A typical diagram for the 77 — >■ AIM process. 



form 



d4>{x,Q) ^ asiQ^) 
ainQ2 



1 



dyV{x,y)cl){y,Q) 



(3.6) 



where V can be computed from a single gluon exchange kernel [46]. While this specifies the 
evolution of (f), the form it has at the starting scale Qq depends on the (non-perturbative) 
details of hadronic binding and cannot be predicted in perturbation theory. However, the 
overall normalisation is set by the meson decay constant /m via [47] 



It was also shown in [46] that for very large the meson wavefunction evolves towards 
the asymptotic form 

^m{x,Q) V^fMx{l-x). (3.8) 

However the logarithmic evolution of (3.6) is very slow and at realistic values the form of 
(j)M can in general be quite different. Indeed, the recent BABAR data [48,49], for example, 
strongly suggests that (l){x,Q) does not have the asymptotic form out to < 40GeV^. 
The issue of the correct form to take for 0(x) has been the subject of much theoretical 
interest (see for example [50]- [54] and references therein), but we will not discuss this in 
any detail here. Instead, to calculate 'benchmark' numerical results we will take 4>{x, Q) 
to have the form proposed in [53] 



It is shown in for example [54] that the available 77 — )■ MM data are quite well described 
by this choice, while this is in general not the case for the asymptotic form (3.8). For 
simplicity we do not include the evolution of (3.9) in our calculation:^ for the realistic 
transverse momenta, k^, values considered in this paper, the scale Q ~ min(x, 1 — x)k±_ 

^It is well known that the annihilation of two energetic (colour charged) gluons into a colourless system 
(two mesons in this case) will in general be accompanied by an intensive bremsstrahlung of relatively soft 
gluons. An exclusive cross section is suppressed by the small probability of not having such radiation, 
which is given by the Sudakov form factor, T, incorporated in the unintegrated gluon densities, /g, in 
(2.1)- see [22, 38] for more details. However, there may also be analogous radiation from the quarks in 
the final state, when the small (~ 1/Q) size qq dipole forms a pion of the 'normal' (~ 0.6 fm) size. The 
corresponding form factor should be accounted for in the evolution (3.6) of the meson wavefunction and in 




(3.7) 



<tP^{x,QQ) = f>V2,fMx{l-x){2x-lf . 



(3.9) 
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is sufficiently low that the effect of including the evolution of (f>{x) is very small, and well 
within other theoretical uncertainties, coming from in particular the choice of 4){x) at the 
starting scale, Qq. 

Finally, returning to the hard amplitudes (3.2) and (3.3), a few comments are in order. 
Firstly, we can see that the amplitudes are divergent as cos 9 — >■ ±1, that is when the 
incoming photons and quark lines become collinear. For small angle scattering, where 
there is no longer a large momentum transferred between the incoming photons and the 
quark lines, the quark propagators become soft and fixed-order perturbation theory can no 
longer be trusted. As we will see in Section 4, upon the inclusion of finite quark masses, 
niq, the amplitudes are finite at cos 9 = ±1, although clearly LO perturbation theory in 
this region can still not be trusted. In practice a reasonable cut-off is simply imposed on 
I cos 6\ to avoid this non-pcrturbative region of phase space: experimentally, for MM CEP 
this simply corresponds to demanding that the mesons are produced relatively centrally, 
and indeed in reality this is always the case. 

Secondly, we can see in the case of neutral meson (ei = 62), for example 7r°7r°, produc- 
tion that the T+^+ , amplitudes (when the incoming photons are in a = state along 
the 77 axis) vanish; we will see in Section 3.2 that a similar result holds for the gg — >■ tttt 
process. For charged pion 7r"'"7r~ production both the \Jz\ = 2 and = amplitudes 
contribute, while in the Tr^vr" case only one term in the jj^j = 2 amplitudes docs, and 
so we would expect the tt^tt" cross section to be suppressed, with an explicit calculation, 
using (3.9) for the pion wavefunction, giving (j(7r°7r'^)/(7(7r"'"7r~) « 0.03 for |cos0| < 0.6 
(see also Fig. 3 of [34]). However, it should be noted that recently BELLE [55] have re- 
ported the significantly large value of cj(7r''7r'^)/cj(7r"'"7r~) = 0.32 ib 0.03 it 0.05 in the range 
\/I = 3.1 — 4.1 GeV and | cos 9\ < 0.6. We would argue that the s values being probed may 
be too small to justify the leading-twist pQCD approximation outlined above (see also [54] 
for some discussion of this) in the tt^tt^ case where the formally leading amplitudes are 
strongly suppressed. For example, if we allow for the fact that the qq pair that form the 
pion can have a non-zero qt then the exact cancellation in (3.2) will not in general occur, 
and this could provide some enhancement to cr(7r°7r°)/cr(7r+7r~). 

3.2 gg — )■ MM: Feynman diagram calculation 

The gg — t- MM process can be calculated in analogy to the 77 reaction, that is using the 
formalism of (3.1), but where the gg qqqq parton-level amplitudes are to be evaluated. 
There are now seven independent Feynman diagrams: four 'abelian' diagrams, where the 
photons are simply interchanged with gluons, and three additional 'non-abelian' diagrams, 
where three and four-gluon vertices are present. For ease of calculation we can group these 
into two separate, individually gauge invariant, sets of diagrams, each weighted with the 
same coIoTir factor, as shown in Fig. 3. In particular, we can decompose the Tgg amplitude 

this way included in the value of (j){x, Q). Moreover, if we choose a relatively low scale <^ Stt-k then we 
have to multiply the result by an additional Sudakov form factor, T(Q^, Snn), describing the probability not 
to emit other gluons with momenta from Q to ^/s^^. Thus we may justify the use of the phenomenological 
form (3.9) only in a limited interval of s^^ ~ 10 — 20 GeV^. For a larger pion pair subenergy we may expect 
an additional Sudakov-like suppression. 
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Figure 3: Basic Feynman diagrams for the gg — )• MM process, grouped into individually gauge in- 
variant subsets Ti (upper) and T2 (lower), and with the relevant colour factors shown schematically. 
The inclusion of all permutations of these diagrams is implicit. 



as 

where Ti and T2 represent the corresponding amplitudes for the diagram sets in Fig. 3. 
Ti is given by the kinematic part (i.e. with the colour factor Cp — 1) of the 77 — t- MM 
amplitudes, (3.2)-(3.3), by setting ei = 62 = 1 and a — )• as- As we have ei = 62, the 
J2 = amphtudes for this vanishes (as in the 77 — >• tt^tt^ case), while for the \Jz\ = 2 
amplitudes we find 

^+__^_+_ 1 -1287r2a|(a - 62)a 

^ ^ sxy{l — x){l — y) a? — h'^cos'^9 ' 

where a, 6 are defined in (3.4)-(3.5). An explicit calculation of the T2 helicity amplitudes 
then gives 

T2++ = T^— = , (3.12) 
_ _ 1 -1287r2a|(a - fe^) cos^ Q 

^ ^ sxy{l — x){l — y) a? — h"^ cos^ 9 ' 

giving for the total gg amplitudes 

^99" = = ' (3-14) 

99 99 Nc sxy{l - x){l - y) a'^ - b'^ cos'^ 6 2 y'''' Nc 

(3.15) 
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the gauge invariance of which has been confirmed by explicit calcTilation. 

We have therefore shown that the gg MM amphtude for = gluons vanishes at 
LO for fiavour-nonsinglet mesons, which, recalhng the = selection rule that strongly 
suppresses the CEP of non-J^ = states, will lead to a strong suppression in the MM 
production cross section, see (2.8). However it should be noted that any NNLO corrections^, 
or higher twist effects^ which allow a = contribution may cause the precise value of the 
cross section to be somewhat larger than the leading-order, leading-twist \ Jz\ = 2 estimate, 
although qualitatively the strong suppression will remain. This is in particular true when, 
as we shall discuss below, the | J^l = 2 amplitude is additionally suppressed by the presence 
of a radiation zero in the amplitude. 

We can sec that the | = 2 amplitude (3.15) vanishes for a particular value of cos^ 6. 
This behaviour, which at first sight may appear quite unusual, is in fact not completely 
unexpected: the vanishing of a Born amplitude for the radiation of massless gauge bosons, 
for a certain configuration of the final state particles is a known effect, usually labelled 
a 'radiation zero', see for instance [56-59] and references therein. It results from the 
complete destructive interference of the classical radiation patterns, leading to a vanishing 
of the amplitude, and the general conditions for the existence of these zeros (which often do 
not occur in the physical phase space region) are derived in detail in [60]. The position of 
the zero is determined by an interplay of both the internal (in the present case, colour) and 
space-time (the particle 4-momenta) variables, as can be seen in (3.15), where the position 
of the zero depends on the choice of meson wavefunction, 0(x), through the variables a and 
b, as well as on the QCD colour factors, see Fig. 4. In particular, the zero occurs in (3.15) 
when COS u 2Cf{o) /Nq, where (a) is the average value of a integrated over the meson 
wavefunctions (j){x) , (j){y) . As < a < 1 for all physical values of x,y and 2Cf/Nc < 1 for 
all Nc (while the prefactors in (3.15) are strictly positive), this will always occur in the 
physical region for any SU(N) gauge theory. 

While this effect, which is present in all theories with massless gauge bosons, is ex- 
pected to occur in QCD, it is usually neutralised along with colour by the averaging of 
hadronisation. The CEP process, for which the fusing gluons are selected to be in a colour 
singlet state by the exclusivity of the event, is therefore in principle uniquely positioned 
to observe these zeros. However, it should again be noted that, as discussed above, as the 
\Jz\ = 2 amplitude is strongly suppressed by the Jz = selection rule, any NNLO or higher 
twist effects which allow a = component to the cross section may give comparable 
contributions; it is therefore not clear that such a zero would in this case be seen clearly 
in the data (at lower values of MM invariant mass, where a non-perturbative framework 
is applicable, we would also not expect to find such a zero, see Section 6). On the other 
hand, the destructive interference effects which lead to the zero in the \Jz\ = 2 amplitude 
(3.15) will tend to suppress the CEP rate. 

'^As there is no contribution from the interference between the 0(a|) 1-loop and the vanishing 0(a|) 
tree level Jz = amplitudes, there arc no NLO corrections in CEP, and the first non-zero perturbative 
correction to the cross section must enter at NNLO. 

*In particular, corrections that result from allowing for a non-zero quark qt with respect to the meson 
momentum, expected to be roughly of order the constituent quark mass qt ~ 350 MeV, may be important 
for the lower values of Mx that are experimentally relevant. 
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Figure 4: Differential cross section dcr/d|cos0| at — 5 GeV, for tlic gg — > AIAI process 
for non-favour singlet scalar mesons. For comparison, the distribution for three choices of meson 
wavefunction are shown, the asymptotic form (j){x) oc x(l — x), the form (3.9) proposed in [53], and 
a (5-function 0(a;) oc 6{x — 5)- 

3.3 gg MM: MHV calculation 

It is well known (see for example [61]) that the tree level n-gluon scattering amplitudes in 
which the maximal number (n— 2) of gluons have the same helicity, the so-called 'maximally 
helicity violating' (MHV), or 'Parke- Taylor', amplitudes, are given by remarkably simple 
formulae [62,63]. These results were extended using supersymmetric Ward identities to 
include amplitudes with one and two quark-antiquark pairs in [61,64,65], where 'MHV 
refers to the case where (n— 2) partons have the same helicity. In these cases, simple analytic 
expressions can again be written down for the MHV amplitudes, while for greater than 2 
fermion-antifermion pairs (recalling that the helicities of a connected fermion-antifermion 
pair must be opposite) no MHV amplitudes exist. More recently, it has been shown in [66] 
that the n-gluon scattering amplitude for any helicity configuration can be calculated with 
this formalism; in particular they can be constructed from tree graphs in which the vertices 
are the usual tree-level MHV scattering amplitudes continued off-shell in a specific way. 
This was extended to include fermions in [67]. 

The basis of all of these results lies in the spinor helicity formalism (see [61,68] for re- 
views): one calculates the matrix elements with the external states having a given assigned 
helicity as an expression written in terms of spinor products of the external particle mo- 
menta. These amplitudes can then be evaluated, usually numerically, squared and summed 
incoherently to give the full cross section. Some basic formulae are given in Appendix A. 

As a further check of (3.14) and (3.15) we can therefore calculate the gg — t- qqqq 
helicity amplitudes using the spinor helicity formalism described above. For the Jz = 

case, the amplitudes (T+_|_ and T ) are MHV and, as we shall show, the vanishing of these 

amplitudes follows simply from the known Parke- Taylor amplitudes. We will therefore 
consider these first, as it is simpler than the \Jz\ = 2 case, and will serve as a clearer 
example. 
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3.3.1 Jz = amplitudes 

In general it is well known that the full n-parton amplitude A4„ can be written in the 
form of a 'dual expansion', as a sum of products of colour factors r„ and purely kinematic 
partial amplitudes 

M„{{pi,hi,Ci}) = ^r„({c^(i)})^„({A;^(i),^^(^)}) , (3.16) 

a 

where {0^(^)1 are colour labels and {^^(i), ^^(i)} the momenta and helicities, respectively, 
of the external legs i = 1 ■ ■ ■ n, and the sum is over appropriate simultaneous non-cyclic 
permutations a of colour labels and kinematics variables. The colour factors T„ arc easy 
to determine, while the purely kinematic part of the amplitude An is to be calculated, and 
encodes all the non-trivial information about the full amplitude, Mn, see for instance [67] 
for more details. 

Wc are therefore interested in calculating the kinematic amplitudes for the 6-parton 

g{±)g{±) — > qqqq process: once this is done we simply use (3.16) to determine the full 
amplitudes. These MHV amplitudes have very simple forms (given in full in [69] and 
elsewhere), and in fact the total n-point amplitude for qqqq plus (n — 4) positive helicity 
gluons can be written down in two lines [61] 



1 {ks ks) {kr k-) 



(kgai)--- {ai kr) {kr h) ■ ■ ■ {hii kg) 

(A"^---A"0im(A''^---A''Oi..2- (3.17) 



Nc {ks ai) ■■■ {ai k-j) {kr h) ■ ■ ■ {hv k-) 

Here the indices r(r) and s(s) refer to the quarks (antiquarks) with colour indices 
and i2(j2)) respectively, and the labels Oj, hi refer to the gluons, while the standard spinor 
contraction '{k,l)^ is defined in (A.l). The sum is over all the partitions of the gluons 
{I + I' = n — 4, I = 0, ■ ■ ■ , n — 4) and over the permutations of the gluon indices, with the 
product of zero A matrices becoming a Kronecker delta and the kincmatical factors equal 
to one when / = (0, n — 4). The overall factor Aq depends on the particular quark helicity 
configuration: for our calculation, it is given by the expressions from [61] for the two quark 
helicity combinations relevant to the meson spin projections. For the equivalent 'MHV 
diagram with all negative helicity gluons, we simply replace {kl) — > [kl], see (A. 2). 

Considering now the case of the 6-parton amplitude relevant to our calculation, we 
make the following identifications 

kr = xk-s kr = {1 - y)ki ks = yk^ ks = {1- x)k3 , (3.18) 
k = 32 12 = ji , (3.19) 

for the collinear quarks (neglecting as usual the qt of the quarks relative to the meson 
momenta) to form colour singlet mesons, where /c3^4 are defined as in Fig. 2. Note that 
there is in general a second possible assignment corresponding to the diagrams for which the 
rr and ss pairs belong to the same mesons, but this does not contribute for non-isosinglet 
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states; we shall discuss this further in Section 4. 

Immediately we can see that first term in (3.17) goes like ~ A;| = 0, while the colour 
factors for the individual pieces contributing to the second term are universal and are given 
by Tr(A"A'') = (5"^/2. Factoring this out, we readily find that the amplitude is given by 



M oc 



+ 



{ki ki){ki k3){k3 k2){k2 k^) {ks ki){ki k2){k2 k^) 
oc {k^ k2){ki ki) + {ki k^){k2 k^) - (fcs ki){ki ^2) = 



+ 



(fcs k2){k2 ki){ki ki 



(3.20) 



from the Schouten identity (A. 4), while the MHV amplitude similarly vanishes. This result 
depends crucially on the colour structure and collinearity of the (massless) quarks/antiquarks 
given in (3.18), which lead to the factorisation of the colour factors and the cancellation 
between the kinematic pieces in (3.20), respectively. It also requires that the produced 
mesons are flavour non-singlet states, see Section 4 for a discussion of this. 

We have therefore in a few lines of algebra confirmed the vanishing of the gg — )• 
MM amplitudes for = initial-state gluons (3.14), which resulted from a non-trivial 
calculation of 7 independent Feynman diagrams. This gives some idea of the power of the 
MHV formalism, which we now apply to the more complicated non-MHV \Jz\ = 2 case. 

3.3.2 \Jz = 2\ amplitudes 




Figure 5: Representative tree diagrams contributing to the g(fci)5(/i;2) — > gggq process with |Jj| = 2 
incoming gluons. Quark labels follow the same notation as (3.17) and ± signs represent particle 
helicity, with all momenta defined as incoming. All contributing amplitudes are of these three types. 



To calculate these amplitudes, which are not MHV, we follow the formalism described 
in [67,70]. The basic idea is that these 'nMHV diagrams can be calculated by connecting 
two MHV diagrams in all allowed ways with a scalar propagator l/p^. For the general 
qqqq case there are four types of diagram that contribute (shown in Fig. 2 of [67]). In fact 
the set of diagrams with a purely gluonic MHV sub-graph can readily be shown to vanish 
when the colour singlet projection is performed and all permutations are summed over (in 
particular, the amplitudes for the diagrams where the external gluon legs are interchanged 
have a relative minus sign but are otherwise identical), and so we will not consider these 
further. Representative diagrams for the three contributing sets are shown in Fig. 5. The 
remaining diagrams are then given by including all possible positions and permutations of 
the two gluons, with the important requirement that gluons are always emitted from the 
same side of the connected quark-antiquark line (see in particular Fig. 1 of [70] and the 
discussion in the text). 
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As described above, this calculation method relies on the factorisation of the total 
amplitude Ain into the sum of products of kinematic partial amplitudes An and colour 
factors Tn- The rules for calculating the kinematic amplitudes are explained in [67], while 
the colour factors are given in [61,67]. For the specific case of the N = 6 parton gg qqqq 
amplitude, the colour factors are given by 

Tn = ^^(A"' • • • A"Onai(A'^ • • • X'ni2a2 (3-21) 

where ii,i2 are the colour indices of the quarks, ai,a2 are the colour indices of the anti- 
quarks and the labels Oj, 6j refer to the gluons. The pair {a} = (0:102) is a permutation of 
the pair {j} = {jij2), where quark is connected by a fermion line to antiquark j^. p is 
then the number of times = jk, with p = 1 if {a} = {j}. As in (3.17), the sum is over 
all the partitions of the gluons (Z + = n — 4, Z = 0, • • • , n — 4) and over the permutations of 
the gluon indices, with the product of zero A matrices becoming a Kronecker delta. For a 
given Tn, the corresponding diagram has the same cyclic ordering of the quark and gluons 
as their colour labels in r„. 

For the two gluon case we can have p = 0, 1 and therefore expect two separate sets of 
diagrams with different colour coefficients to contribute, as can be seen in (3.17) for the 
Jz = amplitude. In particular, recalling the colour singlet assignment of (3.18), the only 
non-zero colour factors are of the form 



2 

The diagrams shown in Fig. 5 correspond to p = 1: the p = diagrams are given by 
making the replacement r s, with the important requirement that a qq pair of the same 
flavour must be connected by a fermion line. We can decompose the previous result for 
the \Jz\ =2 amplitude (3.15) in a similar way 

T^9~ = ^ (^cia - cos^ d) - ^a) . (3.24) 

Following a fairly lengthy explicit calculation, we have then showed numerically that the 
sum of the p = 0,1 amplitudes reproduce these two colour-decomposed terms, thus con- 
firming the result of (3.15). 

We end this section with some final words of explanation. Firstly, to the diagrams 
shown in Fig. 5 we must also add the diagrams where the quarks have positive helicity 
(the helicity of the antiquark must then be negative), corresponding to the second term in 
the meson spin-0 projection; these are obtained by simply interchanging x y. Secondly, 
in implementing the off-shell prescription required to join the MHV sub-graphs, we must 
choose a reference spinor with which to define certain spinor products, see [67] for more de- 
tails. In our calculation we take this to be defined by the gluon momentum \pi). While each 
individual amplitude depends on this choice, of course when the amplitudes are summed 
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this dependence must cancel. However, for certain diagrams unphysical infinities (that is, 
unrelated to standard IR soft/collinear divergences) can arise when external momenta are 
collinear to the reference momentum which, while cancelling in the total amplitude, must 
be dealt with carefully (see, for example [71]). This is achieved by making the replacement 
\Pi) ~^ \Pi) + k)) bringing the divergent amplitudes to a common denominator and then 
using Schouten's identity (A. 4) to simplify the resulting expression, before setting |e) to 
zero, at which stage the result will be finite. Finally, we note that to correctly calculate 
the gg — > MM amplitudes for nonsinglet mesons, care must be taken to omit the set of 
diagrams which only contribute for flavour singlet states, see Section 4. These correspond 
to interchanging the quark legs in the Feynman diagrams in Fig. 3, so that the collinear 
qq pairs forming the mesons are connected by a fermion line, see Fig. 6. We can therefore 
omit these diagrams by requiring that the two qq pairs in the gg — t- qqqq process are dis- 
tinguishable, and using the relevant MHV rules for this case (these are given in [69], for 
example). Using the MHV rules for identical fermions (which includes the permutation 
between the fermion pairs) will not give the required result, as it will implicitly include the 
flavour singlet contribution. 



3.4 gg VV amplitudes 

The amplitudes Tf^^,^ ^^^^ for the g{\i)g{\2) ^(^^3)^(^4) process, where V^(V^) are 
helicity ±1 spin-1 mesons can readily be calculated using the formalism of Sections 3.1 and 
3.2, and are for completeness given here. As we will see in Section 4, for vector meson 
production only diagrams of the type shown in Fig. 3 contribute: in this case, for the 
helicity-0 state, the amplitudes are identical to those for scalar mesons - see (3.14) and 
(3.15). We find 

Tf+^+_ = rf+__+ = rfl^+_ = t5_,_+ = . (3.25) 

n'-'.- = = -w - ^r^'ul i '^'O' - f «) "fir"!' . (3.26) 

^ ^' ^ Nc sxy{l - x){l - y)\ 2 / a^-h^cos^9 

Ncsxy{l-x){l-y)\C'^^ 2 ) a? - cos^ 9 ' 

where it is clear (in the limit that the quark qt = 0) that T^^^^ _|__|_ = Tf^lx^ = ^ 

helicity must be conserved along the fermion line for massless quarks. This result imme- 
diately follows from the helicity conserving gluon-qg vertices which enter the perturbative 
calculation, and the fact that the meson helicity is given by the sum of the helicities of 
its valence quarks: this forms the basis of the so-called 'hadronic helicity conservation' 
selection rule, see [72]. The vanishing of the J2 = amplitudes (3.25) also occurs in the 
77 — >■ VV amplitude, and follows from (3.17) in exactly the same way as for the gg MM 
case^. The | Jj| = 2 amplitudes could in principle be derived using the MHV formalism, 
although we do not consider that here. 



^In fact, in the case of the gg — >■ VV ampUtude, the prefactor Aq in (3.25) also vanishes. 
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4. SU(3)i? singlet contribution: gg ^ riri^r]'r]' ,r}r}' 




In general, as well as the diagrams shown in Fig. 3, it is also possible for the qqqq 
final state to form collinear colour singlet meson states via the process shown in Fig. 6 (a). 
The contribution from this individually gauge invariant subset of 'ladder' diagrams could 
in principle spoil the Jz = cancellation of (3.14). However, for the production of isovector 
states such as pions these diagrams violate isospin conservation and will therefore vanish: 
for example, for -k^-k^ production the \uu) and \dd) contributions interfere destructively to 
give zero total contribution. 

Calculating the relevant amplitudes in the usual way we find for scalar mesons, for the 
valence quark contribution (these results could also be derived using the MHV formalism 
described previously, although we do not consider it here) 



647r2a| (l + cos^^ 
Nc sxy{l — x){l — y) (1 — cos^ ^^"^ 



^lad. ^ ^lad. ^ ^ 647r2a| (l + 3cos2g) 

+- -+ Nc sxy{l -x){l- y) 2(1 - cos2 Of' ^ ' ' 

While these amplitudes do not contribute, for example, for vrvr production, they will in 
general be relevant for r/'r/' production and, through rj — rj' mixing, rjrj production. The 
SU(3)i? basis states are defined as 

\r]o) = — -=|uu + dd + ss) , (4-3) 
v3 

\t]s) = -^\uu + dd - 2ss) . (4.4) 

In terms of these states the rj and rj' wavefunctions are given by 

I?;) = cos6'p|??8) - sin6'p|r/o) , (4.5) 

\r]') = sm9p\ris) +cos9p\r]o) . (4.6) 
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In the case of rji] CEP, the flavour singlet state |ryo) receives a non-zero ,7^ = contribu- 
tion through (4.1), and may in general greatly enhance the cross section relative to tttt 
production. There is some variation in the values of 6p extracted from data, depending 
on the theoretical assumptions used in the analysis: following chiral perturbation the- 
ory, for example, the meson masses in the SU{3)f nonet can be used to extract a value 
9p = 11.5° [73], but from fits to the decays of different light vector and pseudoscalar mesons 
(see for instance [74,75] and references therein), a value of 9p ~ —15° is favoured, and we 
will take this value as a default^*^. In this case, the flavour singlet contribution will be 
suppressed by a factor sin^ ~ 1/200, which may therefore be comparable to the | J^l =2 
flavour-octet contribution. Moreover, the flavour-singlet amplitude is further enhanced by 
the larger multiplicity of flavour states that can contribute to the independent qq pairs 
in Fig. 6: for each meson, there are 3 contributing flavours (mm, dd and ss), enhancing 
the amplitude by a factor of 9, while the wavefunction normalisation in (4.3) only gives a 
factor of 1/3. The amplitude squared is therefore enhanced by a factor of 3^ and we may 
therefore expect, in regions of phase space where the perturbative formalism is applicable, 
the 7777 CEP cross section, despite the mixing angle suppression, to be dominant over tttt 
CEP - we shall see in Section 5 that this is indeed the case. 

As the 1]' state is dominantly flavour singlet, we also predict that rj'-q' CEP should be 
strongly enhanced relative to vrvr and 7717 production. We note that the gg — t- rjrj' process, 
which vanishes for the diagrams given in Fig. 3 due to the orthogonality of the r] and r/' 
flavour wavefunctions, can also occur via the diagrams of Fig. 6 (a), although the cross 
section is predicted to be suppressed by a factor of sin^ 9p relative to rf'rf' production. 

It also immediately follows from helicity conservation along the quark lines that the 
production of transversely polarized vector mesons cannot proceed via these diagrams (in 
the limit that the quark qt = 0). Moreover, for longitudinally polarised vector mesons, we 
find 

T'^tfiO = = (2a - 1) T|t , (4.7) 

r|-f ,00 = T'^tfiO = (1 - 2a) r|-f • , (4.8) 

which are antisymmetric under the interchange x o (1 — x) (or y -H- (1 — y)), and will 
therefore vanish upon integration over the (symmetric) meson wavefunction, (/>(a;). This 
result for these 'ladder' diagrams, where each meson state couples separately to two gluons, 
recalls the well-known Landau- Yang theorem [76], which states that a spin-1 particle cannot 
couple to two on-shell massless vector bosons^^: although the intermediate t-channel gluon 
is in general far off-shell, the overall amplitude nevertheless vanishes. It therefore follows 
from (3.25) and (4.7)-(4.8) that the CEP of light vector pairs will be strongly suppressed 
in the perturbative regime, irrespective of their flavour structure, and indeed the pp, lou 
and (f)(f) rates are predicted within this formalism to be the same up to small (higher order, 

'^'^Often in the literature, rather than using the singlet-octet basis given in (4.3) and (4.4), a 'quark' basis 
is taken, with corresponding mixing angle 4>. This is related to 9p via Op = (j) — arctan v^. 

^^In the hard exclusive formalism, the V —¥ gg amplitude is proportional to (2a; — 1), and will therefore 
vanish upon integration over (f){x). 
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higher twist) corrections. 

Finally, wc should in general consider the possibility of a two-gluon Fock component 
\gg) to the rj{ri') mesons (see for example [77,78]). The lowest order {a^g) diagram, where 
the incoming gluons couple directly to the meson states, with an intermediate i-channel 
gluon exchanged between the mesons, will be strongly suppressed by the virtuality of the 
exchanged gluon, that is by the 77(77') gg* form factor, where the g* is far off-mass-shell. 
In the hard exclusive formalism outlined in the previous sections, these form factors are 
modeled perturbatively: a representative diagram of the valence gluon contribution to the 
gg —7- ri{r]')r]{ri') process is shown in Fig. 6 (b). Thus it occurs to the same order in as 
as the valence quark process, with the size of the gluon contribution determined by the 
size of the meson two-gluon wavefunction 0g(a;). Unfortunately, the size of the two-gluon 
component of the 77(77') wavefunction is very poorly determined from experiment, and there 
exists no firm theoretical consensus as to its relative importance^^. While as — >■ 00, it 
can be shown that the gg wavefunction vanishes due to QCD evolution [79, 80] 

hm 0g(x) = O, (4.9) 

there is no reason to assume this will be the case at the experimentally relevant energies. 
In [78], to consider one example, a fit to the r]{r]')^ transition form factor, F^(^/)(Q^), 
suggests that the size of the gg wavefunction may be of the same order as the quark con- 
tribution, although the extracted value is consistent with zero. Moreover, this fit assumes 
that the qq wavefunction is close to the asymptotic form, which appears to contradict the 
more recent BELLE data [48,49], see also [54]. There is also some indication from Xc 
decays that the role/contribution of the '55' component may be smaller than that found 
in [78]. In particular, no enhancement in the XcO ~^ v'v' relative to the XcO 7/r/ branching 
is observed, while the Xc2 — ^ v'v' decay is in fact strongly suppressed [81]: this cfi'ect would 
be quite surprising if a sizeable 'yy' component to the 77' wavefunction were present. While 
various other estimates are also available in the literature (see for example, [74,77,82-84]) 
no firm consensus exists about the precise size of the gg contribution. 

Without a reliable input for (t)g{x), we therefore do not consider a full numerical cal- 
culation of the gg contribution here, although this can readily be performed using the 
formalism outlined in the preceding sections. As described above, the relevant perturba- 
tive gg — t- 4g' amplitude can be calculated in the usual way^^ and, as this does not vanish 
for J2 = initial-state gluons, we may expect it to enhance the 7777 and rj'rj' CEP rates 
(although, as with the \qq) flavour singlet contribution, we will expect the gg contribution 
to the [77) state to be fairly small). This will depend sensitively on the size of the two-gluon 
wavefunction: therefore, by considering the CEP of 77(77') pairs at sufficiently high invariant 
mass, it may be possible to extract some information about the relative importance of the 
leading-twist quark and gluon wavefunctions. 



We are very grateful to V. Chernyak and A. Groznin for a discussion of these issues. 
The relevant gg — >■ ggqq amplitudes should in general be included as well. 
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5. MM CEP: results 

To calculate the MM production cross sections, we use the formalism described in Sec- 
tion 2. We take /jr = 133 MeV and = fp = 200 MeV, as in [35], and assume the flavour 
octet and singlet decays constants for states (4.3) and (4.4) are given by fo = fs = fn- 
although this will in general not be true [75], we consider any deviation from this (which 
occurs at the 10 — 20% level) to be within the uncertainties of the calculation. We assume 
in all cases a universal meson wavefunction given by (3.9). 

We first show in Fig. 7 the spin-averaged differential cross section dfi/dl cos(0)| for the 
gg — )• MM subprocess (the distributions are of course symmetric under the interchange 
cos — )• — cos 9). In all cases the cross sections are strongly enhanced as | cos 0| — )• 1, where 
they are in fact divergent in the limit of massless quarks taken here. Clearly, in the limit of 
small angle scattering we can no longer trust the calculation, but out to reasonable values 
of cos 9 we expect the dominant angular behaviour to be described by these distributions: 
indeed measurements at BELLE of 77 — t- vr+vr^ (and 77 — )■ K^K^, which is expected 
theoretically to have a very similar angular distribution) out to |cos0| = 0.6 [81] and 
77 — )■ vr'^vr'^ out to |cos0| = 0.8 [55] observe the 1/sin^^ distribution predicted within 
the perturbative framework. It should also be noted that, upon the inclusion of non-zero 
quark masses, the amplitudes are not in fact divergent as |cos0| — t- 1. An interesting 
example of this is in the case of the | = 2 flavour singlet amplitude (4.2), which, due 
to conservation of the projection of the total angular momentum, J, onto the z-axis, must 
vanish in the limit of zero angle scattering. While for massless quarks this is lost due to 
the IR divergence of the quark propagators, upon the inclusion of a non-zero quark mass 
the amplitude does indeed vanish at | cos9\ = 1, see Fig. 8. Although we do not show it 
explicitly here, a similar result must of course hold for the nonsinglet amplitude (3.15), and 
we note that the non-zero gg — )• VV helicity amplitudes (3.26, 3.27) vanish as expected 
in the forward or backward limit, when the initial and final-state spin projections on the 
2;-axis are opposite. Nonetheless, we of course cannot trust fixed-order perturbation theory 
in this region of phase space. 
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Figure 7: The differential cross section da/dj cos(6')| for the spin averaged gg — MM process for 
various meson states. 
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Figure 8: gg MM amplitude for production of flavour singlet mesons (4.3), with gluons in a 
\Jz\ = 2 state. Plotted is the amplitude as a function of | cos(6')| for massless quarks, and with the 
inclusion of a non-zero quark mass. For the purpose of illustration we choose the value of vriq = 0.35 
GeV, but it should be noted that the amplitude vanishes at | cos 6*1 = 1 for any non-zero quark mass. 



It should be noted however that the angular distributions in Fig. 7, which are averaged 
over all incoming gluon helicities at the cross section level, are not directly relevant to the 
case of meson pair CEP, and are only shown for illustration. In particular, we have noted 
the importance of the = selection rule that operates for CEP. This has significant 
consequences for the case of MM CEP: in particular, the subprocess amplitudes for the 
flavour- nonsinglet scalar (3.15) and all vector (3.25) states vanish for = gluons, and 
we therefore expect the CEP cross sections to be heavily suppressed, see (2.8), relative to 
the cross section for flavour-singlet states {rj'rj' , rjT], rjr]'), where the subprocess amplitudes 
do not vanish for = incoming gluons, see (4.1). The r]r] and r]r]' CEP cross sections 
are strongly dependent on the precise level of rj — ij' mixing, through which a = 
component can enter. In Fig. 9 we show the CEP cross sections da/dMx, where Mx is 
the invariant mass of the meson pair, for the production of various scalar and vector states. 
The suppression of the tt'^tt'^ and vector meson cross sections is clear^'*, in particular in the 
TTTT case where the radiative zero in the angular distribution will tend to further reduce 
the cross section. The r]'r]' cross section, on the other hand, is predicted to be quite large. 
The vector meson pp cross section is also shown: within the perturbative formalism, the 
and COLO rate are to lower order expected to be identical, see Section 3.4. 

In Fig. 9 we also show the CEP cross sections for the production of the saixie states, as a 
function of the cut, Ecnt, on the meson transverse energy E± = E sin 9 (where 9 is the angle 
of the pion from the beam axis in the lab frame), and the same effects are clear. In Fig. 10, 
we show the effect of increasing the c.m.s. energy, ^/s, and cut on the meson pseudorapidity, 
r/jvf ) on the cross section, taking the case of rf'-q' CEP for illustration. Unfortunately, for 
the low X and values probed in the CEP of lighter mass objects, there is a large degree 
of uncertainty in the single PDFs. Recalling the strong PDF dependence (ctcep ~ {xg)'^) 
of the CEP cross section, the effect of this will be quite severe. This is already relevant 

^*The charged tt^tt" and p~ CEP cross sections are expected to be a factor of 2 larger from isospin 
symmetry and the non-identity of the final state particles. 



- 21 - 



at Tevatron energies, but is a stronger source of uncertainty when it comes to making 
predictions for the LHC, where the probed x values are even lower. In Fig. 11 we show this 
explicitly by comparing the predicted vr^vr'^ CEP cross section as a function of -Ecut for two 
different choices of PDF set, MSTW08LO [85] and MRST99 [86], as in [8] for the case of 77 
CEP. Clearly this is a significant source of uncertainty in the normalisation of the predicted 
cross sections, with the MSTW08LO and MRST99NLO sets providing approximate upper 
and lower bounds, respectively, on the range of predictions coming from different PDF sets. 
We use these two particular sets because we would argue that they span a realistic range 
of small X parton distributions. MSTW08LO is the outcome of a leading-order pQCD 
global fit to an up-to-date set of DIS and other hard scattering data. At small x, it gives 
a reasonable, though not perfect, description of HERA F2{x, Q^) data. The corresponding 
NLO version (MSTW08NLO) has a very different gluon distribution at small x and Q^, 
with g{x, Qq) < for X < 10~^. Because of this behaviour, the MSTW08NLO gluon gives 
unstable results when used to calculate the skewed PDF, fg. We prefer to use the older 
MRST99 NLO set, which has a more benign small-x form, while still retaining the essential 
features of a NLO fit, in particular a gluon that is smaller at small x than at LO. We take 
the MRST99 set, which gives predictions that are in better agreement with the existing 
Tevatron data, as our default set, see [8] for more discussion of these issues. 
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Figure 9: da /din M'^ for meson transverse energy E± > 2 GeV, and cross section as a function 
of the cut Ecut on the meson E± at -y/s = 1.96 TeV for the CEP of meson pairs, calculated within 
the perturbative framework, including both the Jz = and | J^j =2 contributions. The mesons are 
restricted to have pseudorapidity |77a/| < 1.8. Meson masses are neglected throughout. 



As was discussed in the Introduction, the calculation of the Tr'^vr^ CEP cross section 
has important consequences for the possible vr'^vr'^ background to 77 CEP, when one photon 
from each tt^ decay is undetected or the two photons merge. At first sight it would appear 
that the cross section for this purely QCD process could be much larger than that in the 
77 channel and so would constitute an appreciable background, but fortunately this is 
not the case. Firstly, we have seen that the amplitude to form an exclusive pion with 
large transverse momentum, k±, is proportional to the ratio fn/\^ ~ fn/k± (see (3.7) and 
(3.15)), that is the cross section of the gg — t- vr'^vr'^ hard subprocess contains the numerically 
small factor {f-w/k^)^ which in the region of interest is comparable (or even smaller) to 
the QED suppression, {o-qed / 013)"^ , of the gg — )■ 77 cross section. Secondly, and crucially. 
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Figure 10: rj'ri' CEP cross section at y/s — 1.96 TeV, as a function of tlie cut E^ut on the 
meson E±, calculated within the perturbative framework, including both the Jz = and \Jz\ —2 
contributions. Plotted is the cross section for different values of c.m.s energy y/s and cut on the 
meson pseudorapidity r]M- Meson masses are neglected throughout. 
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Figure 11: tt^tt" pair CEP cross section as a function of the cut Ecut on the meson E±, at 
\/s = 1.96 and ^/s — 7 TeV, calculated within the perturbative framework, including both the 
= and I J^l = 2 contributions, for two choices of PDF, MSTW08LO and MRST99. 



the vanishing of the LO amplitude gg — )■ vr vr with = initial-state gluons leads to 
a further ~ two orders of magnitude suppression in the CEP cross section. We therefore 
expect the background contribution to 77 CEP to be smah ° - we shall address this 
issue in more detail in a forthcoming publication [87]. 

Our results can readily be extended to the kaon sector, although we do not consider this 
numerically here. In particular, under the assumption of exact SU (3) flavour symmetry, 
the K^W and K+R- CEP cross sections can be calculated in the same way as the vr'^vr'^ 
and vr^vr^ cross sections, with the replacement /tt — ?• fx- In fact, SU{3) flavour symmetry 
breaking effects can be non-negligible, and to precisely estimate the K^K and K~^K~ 
cross sections, a modified narrower form of the meson wavefunction, which accounts for 
asymmetry between the s and (u, d) quark masses, should be taken, see [35] for more 
details. Without the inclusion of this modified wavefunction, the perturbative formalism 

'^^As noted in Section 3.2, the cross section may however be somewhat enhanced by higher-twist or 
NNLO corrections which allow a = component. At lower values of tttt invariant mass, an additional 
'non-perturbative' contribution must also be considered, see Section 6. 
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tends to overestimate the 77 — >• K^K~ cross section when compared to BELLE data [81]. 

FinaUy, we note that the CEP of these states is now included in the pubhcly available 
SuperCHIC Monte Carlo [88], which also generates the CEP of a range of low-mass states 
(Xc,b, rich, 11) ■ 

6. Non-perturbative meson pair production 




Figure 12: Non-perturbative meson pair [M^, M4) CEP mechanism, where M* is an intermediate 
off-shell meson of type M . 



For low values of the meson pair invariant mass we may not expect the perturbative 
framework described in Sections 3.1 and 3.2 to be applicable. In this region of phase space 
we can model the process using a 'non-perturbative' picture [30, 32, 89] shown in Fig. 12: 
the MM pair are created via double-Pomeron exchange, with an intermediate t-channel 
off-shell meson, M*. 

We can calculate the resultant amplitude using the standard tools of Regge theory, see 
for example [90]. The CEP cross section is given by 

= I dpl,dpl,dy^dy,dkl\-^ , (6.1) 

where ^/s is the c.m.s. energy, pij_,p2i. are transverse momenta of the outgoing protons 
and k± is the meson transverse momentum. is the soft survival factor, expected to be 
of a similar size to Sgj^, in the perturbative case (2.6). The matrix element is given by 
Ai = Aif + Aiu, with t = {Pi — ^3)^, u = {Pi — A;4)^, where Pi is the momentum transfer 
through Pomeron i, and ^3^4 are the meson momenta. We have 

^^=Jih'MJFM.)n.4{^) (If) , (0.2) 

where M is the meson mass and we take sq = 1 GeV^ and ap(p?_|_) = 1-08 — 0.25]p?_|_], for 
measured in GeV^ [91], and Sij = {p[ + kj)"^ is the c.m.s. energy squared of the final 
state proton- meson system {ij). The proton form factors are as usual taken to have an 
exponential form, Fp{ti) = exp(5jtj/2), while the Pomeron trajectory is included in the 
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definition of the slope 

Si = 60 + 2a' log (^^^ , (6.3) 

with 60 = 4 GeV~^. Concentrating on the case production, the overall cross section 

normalisation is set by the total pion-proton cross section a{pTr) = (Joisij / so)°'^^'>~^ 30 
mb at the relevant sub-energy. 

Finally, we have to include an additional suppression factor to calculate the genuinely 
exclusive MM cross section - the small probability not to produce any additional hadrons 
in the Pomeron-Pomeron fusion process, which can be evaluated by taking a Poisson dis- 
tribution in the multiplicity of secondaries. Note however that we cannot assume a Pois- 
son distribution in the multiplicity of all charged tracks, as this ignores the correlation 
N'^ = N~ . Secondly, a fraction of the secondary pions is produced via resonance decay 
or even in a minijet fragmentation - in the latter case wc should consider a Poisson distri- 
bution in the number of minijets. Finally, besides charged particles there are also neutral 
secondaries to consider. For the above reasons we assume a Poisson distribution in the 
number of negatively charged particles. The mean multiplicity, n = {N~), grows with 
the Pomeron-Pomeron energy, s = {Pi + ^2)^, as n ~ c ■ ln(s/so), with the coefficient-'^^ 
c ~ 0.5 — 1. The Poissonian probability not to create any additional secondaries is then 
given by w = exp(— ri) = {s/sq)~'^. For numerical estimates we take c = 0.7 here. 

One may ask what the origin of this additional suppression, w, is in terms of the non- 
perturbative diagrams of the type shown in Fig. 12. We should recall that this is just a 
representative diagram: besides 'eikonal' proton-proton rescattering (not shown in Fig. 12, 
but included as the soft gap survival factor = S'^^^ in (6.1)) we have to account for the 
possibility of additional meson-proton (and meson-meson) interactions which play the role 
of the survival factor S'^^j^ shown in Fig. 1, and for the reggeization of the virtual t-channel 
meson M*. It is this reggeization of the meson M* trajectory which leads to the additional 
power suppression described above. In the case of tttt CEP for example, with the slope of 
the pion trajectory taken as = 0.9 GeV~^, the coefficient c = 0.7 corresponds to a mean 
momentum transferred —t = 0.4 GeV^. Actually, here we are considering the case of pions 
with a rather large k± > 1 — 2 GeV, for which the contribution to the diagram shown 
in Fig. 12 from pion reggeization, taken literally, would be negligible. However besides 
the pion contribution in the place of M* there may be 02 and/or ai Reggc trajectories, 
and the whole structure of this non-perturbative interaction is therefore very complicated. 
To get an order of magnitude estimate of this effect we consider the simplified diagram of 
Fig. 12, neglecting the pion reggeization but including the probability, w, not to produce 
any additional secondaries, the procedure for calculating which wc have described above. 

We note that the CEP of vrvr pairs within the non-pcrturbativc framework was consid- 
ered recently in [36] (see also [32,37]). The model used is largely similar to that considered 
here, but this extra suppression is not included and so their estimates for the tttt non- 
perturbative cross section in general become significantly larger than our predictions as 
M;r7r is increased, leading for example to roughly an order of magnitude difference at the 

^®It would be interesting to measure the height of the 'plateau', dN/drj, produced in Pomeron-Pomeron 
(DPE) collisions at the LHC. 
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Figure 13: Comparison between perturbative (2.1) and non-perturbative (6.2) tt+tt^ CEP cross 
sections as a function of tlie cut E'cut on tlie pion transverse energy E± . In the non-perturbative 
case, an exponential pion form factor Fjr{t) — exp{bt), with 6=1 GeV^^ is taken. 




I I I I I I I I I J 1 L I I I I I I I 

3.2 3.3 3.4 3.5 3.6 3.7 3.8 ' 3.2 3.3 3.4 3.5 3.6 3.7 3.8 

M„ [GeV] M„ [GeV| 

Figure 14: Comparison between perturbative (2.1) and non-perturbative (6.2) tt+tt" CEP dif- 
ferential cross sections da/dM, for different cuts on the pion transverse momentum, k±, and for 
the pion pseudorapidity restricted to lie in the region i] < 1. In the non-perturbative case, an 
exponential pion form factor FT^{t) = exp(6t), with 6=1 GeV^^ is taken. 



Mtttt ~ mass.^^ Moreover, the authors of [36] do not consider the possible perturba- 
tive contribution due to the gg — t- vrvr subprocess, which we expect to be dominant, and 
therefore increase the overall rate, for higher values of M-^^. In fact, these two differences 
may compensate each other in the M-,^-,^ ~ mass region, leading to an approximate 
consistency between our results. 

In the case of non-perturbative MM production, large meson k±_ values are suppressed 
by the form factor, Fjuit), of the intermediate off-shell meson in (6.2), although in this 
case it is not completely clear what form to take for F/v/(i)- In particular, we may either 
take the 'soft' exponential exp(6t), or the 'hard' power-like form ~ 1/t, both of which are 
used in the literature [31,32,89]. As |t| (that is, the meson k'^) is increased, these two 
choices of form factor give vastly different cross sections, with the power-like form factor 

^^In the previously measured mass range (M^^ < 2 GeV) at the CERN ISR [92], where the effect of this 
suppression is not too large, our results are roughly consistent, when we account for the secondary Reggeon 
contributions included in [32], but which we have omitted above for simplicity, as their contribution will be 
insignificant at the c.m.s. values we consider here. 
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giving a significantly larger rate. However, in general we should expect a smooth matching 
at a reasonably low value of kj_ between the non-perturbative and the perturbative region 
of phase space, where we expect the mechanism outlined in the previous sections to be 
dominant. At higher values of k_i_, the contribution from the power-like form factor should 
not be included, since it represents a form of double-counting of the perturbative tail region, 
which we should calculate within the pQCD CEP formalism. We show this in Fig. 13, 
where the non-perturbative (with exponential form factor) and perturbative 7r"'"7r~ CEP 
cross sections are plotted as a function of the cut on the E± of the pions: the matching 
between the two regions occurs at the reasonable value of E± ~ 1.5 GeV. Due to the 
strong Jz = suppression of the perturbative cross section, the transition region occurs at 
a somewhat higher value of pion E± (and hence Mtttt) than we would otherwise find, and 
this has the consequence that both the perturbative and non-perturbative mechanisms may 
be relevant in the Mj^^ ~ My. mass region relevant for evaluating the potential continuum 
background to x — ^ tttt CEP. Plotted in Fig. 14 is the differetial cross section da/dMT^T^ for 
the perturbative and non-perturbative contributions in the mass region relevant to Xc CEP 
at -y/s =1.96 and |ry^| < 1. When no cut is placed on the pion transverse momentum the 
non-perturbative mechanism is predicted to be dominant in the Xc mass region, however 
when more realistic experimental cuts are placed on the pion k±, the relative contribution 
of the non-perturbative mechanism is supressed by the pion form factor, -F^(t): this can be 
seen for example in the case of the cut k± > 1.5 GeV, for which the perturbative mechanism 
is predicted to be dominant. We note that, when available, CDF data [5, 93] on central 
exclusive vr^vr^ production may provide useful information about the off-shell pion form 
factor and help further illuminate this issue. 

In fact, the CEP mechanism shown in Fig. 1 does not represent the true perturbative 
tail to the double Pomeron exchange process shown in the Fig. 12: we shall consider this 
process in the following section. 



7. MM CEP: secondary mechanism 
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Figure 15: (a) Representative perturbative diagram for 'skewed' meson pair CEP. (b) Represen- 
tative perturbative diagram for 'symmetric' meson pair CEP. 
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As well as the standard CEP diagram shown in Fig. 15 (a), there is in general a 
second mechanism^^ for producing meson pairs at high k_\_. This is shown in Fig 15 (b): 
the second f-channel gluon now couples directly to a quark line, with the collinear qq 
pairs forming mesons in the usual way, in both the flavour-nonsinglet and flavour-singlet 
combinations described in Sections 3 and 4, respectively. We shall label Figs. 15 (a) and 
(b) as 'skewed' and 'symmetric' CEP, respectively: in Fig. 15 (a) the second t-channel 
gluon is much softer than the fusing gluons, whereas in Fig. 15 (b) both t-channel gluons 
participate symmetrically in the hard subprocess. The symmetric diagram represents the 
perturbative tail to the non-perturbative double Pomeron exchange mechanism shown in 
Fig. 12, resolving the two-gluon structure of the exchanged Pomerons. As this diagram is 
of the same order in 0:5, it might in principle give an important contribution to the meson 
pair CEP cross sections given in Section 5. Moreover, in the case of scalar non-singlet 
(and vector meson) skewed CEP amplitude, which are suppressed by the = selection 
rule, the effect may be particularly important: in Fig. 15 (b), the vrvr and vector meson 
production amplitudes will no longer vanish for forward outgoing protons. 

The symmetric CEP amplitude, considering for simplicity the case of forward outgoing 
protons {pi_^ = 0), can be written in the form 

Tsym. =Vr^/ 4 ■M.syra.fg{xuXi,ql^,l?;ti)fg{x2,X2,ql^,l?\t2) , (7.1) 

where the notation follows from Fig. 15 (b) and the subprocess amplitude TWgym. is given 

by 

A^sym. = ^^^'^'^^'^''^i'x52±'??±52\Ca'^ (7-2) 
- ^V^sym. . (7.3) 

For example, considering the specific case of scalar non-singlet meson production, T^ym., as 
shown in Fig 15 (b), is given by 

l/sym. = 167r «5 2]v^ j dxdy(p{x)(p{y) • (7.4) 

The other diagrams corresponding to interchanging the outgoing quark legs l\ O I2 and 
^3 ^4; as well as the s-channel diagrams containing 3- gluon vertices, should also be 
included, and the resulting expression can then be combined with (7.1) to give an explicit 
evaluation of the full symmetric CEP amplitude. 

However, simply by inspecting the form of the amplitude (7.2)-(7.4) we find that the 
symmetric process in general represents a power correction to the skewed CEP process, 
and therefore formally gives a subleading contribution. We demonstrate this in full in 



The authors are grateful to JefT Forshaw for bringing this type of diagram to our attention. 
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Appendix B, where we show that 



^skew. "'_|_ 

where k± is the meson transverse momentum, and {Q'jJ is the average gluon transverse 
momentum in Fig. 15(b), see (B.ll). That is, the symmetric amphtude is suppressed by 

the meson A:^ ~ s, and its relative contribution will decrease as the meson k± is increased. 
The skewed CEP amplitude is therefore, at least formally, dominant at high k±. However, 
if we take k± « Mx/2 = 5/2 GeV we find 

(Tsvn,. 1 1 



^skew. 4 2 

for a typical {Q\) = 3 — 4 GeV^. At smaller values of k±, it therefore appears that the 
suppression of the symmetric contribution may not be numerically too large, although as 
fe_L is increased we will indeed expect the symmetric contribution to the CEP cross section 
to be negligible (falling as ~ This is only a rough estimate which is based on 

kinematics alone, with the precise result requiring a specific calculation, but, when we 
recall that for example the skewed amplitude for tttt production is strongly suppressed by 
the Jz = selection rule (by ~ two orders of magnitude), it appears that the symmetric 
diagram may indeed be important for lower values of Mx, and this is equally true for vector 
meson production. On the other hand, for non-suppressed rj'rj' CEP, and the relatively non- 
suppressed rjT]' and rji] CEP, we can safely neglect the symmetric contribution. 

We note that the skewed amplitude shown in Fig. 15 (a) is almost purely imaginary 
while the non-perturbative amplitude in Fig. 12 and the symmetric amplitude in Fig. 15 
(b) are almost purely real. In the case of the non-perturbative amplitude, there are two 
Pomeron exchanges: each Pomeron has intercept Q!(0) ~ 1 and therefore gives a dominantly 
imaginary contribution, leading to the dominantly real product i'^ = —1. In the case of 
the symmetric perturbative diagram we have two gluon loops {qiqi and q2Q2), with each 
loop representing the high sub-energy, + li + ^4)^ or (j)'2 + ^2 + ^3)^) positive signature 
amplitude, which is equivalent to a Pomeron exchange and gives an imaginary contribution; 
the total amplitude is therefore real.^^ On the other hand, the skewed diagram in Fig. 15 
(a) contains only one gluon loop (with loop momentum Q±), which therefore gives an 
imaginary contribution. 

We stress that the estimate of (7.6) is only approximate: there may therefore be 
important physics that this misses. Firstly, we have ignored the effect that including 
the exact form of the subprocess production amplitudes (see for instance (7.4)) may have 
on the relative size of the skewed and symmetric contributions, although it can readily 
be seen that in terms of, for example, colour structure, the symmetric amplitude is not 
further suppressed. Secondly, and more importantly, is the issue of how to include the 



^®The main features of the symmetric perturbative contribution follow closely those of the non- 
perturbative contribution. The perturbative amplitude of Fig. 15 (b) may therefore be considered as the 
large k± tail of the non-perturbative amplitude. 
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generalised PDFs in (7.1). In the case of the skewed CEP amphtude shown in Fig. 15 
(a), the skewed PDFs can be related to the standard integrated (or unintegrated) PDFs 
(see Section 2) by using the fact that the light cone momentum fraction carried by the left 
screening gluon, x', is much less than that, x, carried by the active (right) gluon, \x'\ <^ x. 
On the other hand, for the symmetric case of Fig. 15 (b) both the exchanged gluons have 
comparable momenta, qi ~ g,. Moreover, in contrast to the diagram for the DIS cross 
section, where the momentum fractions carried by qi and qi gluons have different signs 
(with the directions indicated as in Fig. 15 (b)), here they have the same sign. In terms 
of the variables conventionally used to describe the generalised PDFs, X = {x — x')/2 and 
( = {x + x')/2, this corresponds to the so-called 'time-like' \X\ < ( domain, where the 
skewed parton distribution cannot be extracted from DIS data, which lie in the 'space- 
like' domain. In the |X| < domain the generalised PDFs describe the wavefunction of a 
'gluebair formed by two t-channel gluons. Since we do not know the wave (distribution) 
functions of the appropriate 'glueballs' and their couplings to the proton, the best we can 
do is to put an upper limit on the cross section, based on the Schwarz inequality [94] 

fg{x, x', Q^ ...) < ^{{fg{x, -X, Q\ ...) + fgix', -x', ...)) , (7.7) 

where the diagonal gluon distribution fg{x, —x,Q'^ , ...) can be extracted from DIS data. 
Finally, we note that in the symmetric amplitude we must be careful to forbid perturbative 
emission from both t-channel exchanges, which can both be hard, although in general at a 
scale that is somewhat lower than /x = Mx/2. We leave a full quantitative investigation of 
these issues to a forthcoming publication [87] . 

8. Conclusions 

In this paper we have studied the central exclusive production of meson pairs, MM, for a 
variety of light meson states at hadron colliders. We have concentrated on the CEP process 
in the perturbative regime, that is where the MM invariant mass is sufficiently high that the 
formalism outlined in Section 2 can be applied, with the gg — t- MM subprocess described 
by the 'hard exclusive' formalism outlined in [34,35,46]. In particular, the gg — )■ MM 
helicity amplitudes are calculated for a range of light scalar {it'^'k~, K^K~, K^K , r/r/, 
r/'ry', ryr?') and vector (p/9, ujuj, (fxf)) meson states, and the full MM CEP cross section is 
then calculated using the formalism outlined in Section 2. The CEP of these states has 
also been included in the publicly available SuperCHIC Monte Carlo [88], which generates 
the CEP of a range of low-mass states {xc,b^ Vc.bj ll)- 

We have observed and discussed various interesting theoretical properties of the gg — )• 
MM helicity amplitudes. In Section 3.2, we have seen for example that the LO tttt produc- 
tion amplitude (and more generally that for any pair of flavour non-singlet scalar mesons, 
e.g. K'^K~ , K^K^) vanishes when the fusing gluons are in a = state. This vanishing 
also occurs for the production of vector mesons (e.g. uou, (p(f>), sec Section 3.4. These 
results can be considered as generalisations of the vanishing of the known 77 — ^ MM 
amplitudes for = incoming photons, for neutral scalar and neutral and charged vector 
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mesons, see [46]. We have also seen that the \Jz\ =2 amphtude for the production of scalar 
mesons contains a radiation zero, where the LO amplitude vanishes exactly for a particular 
value of cos^ 9 that depends on the colour factors of the SU (3)c gauge group as well as the 
choice of non-perturbative meson wavefunction, ^{x); this destructive interference effect 
will tend to suppress the production cross section. 

We have also applied the tools of the MHV formalism [61-63] to the gg — )■ MM 
calculation. In particular, the = vanishing of the scalar meson amplitude, which 
follows from a non-trivial calculation involving many independent Fcynman diagrams, has 
been shown to follow simply from the known Parke- Taylor amplitudes, while the | J^:! = 2 
results have been confirmed using a generalisation of MHV techniques to the calculation of 
non-MHV amplitudes, see [66,67]. While in this paper we have only explicitly calculated 
the amplitudes for scalar flavour non-singlet meson production within this formalism, it 
could readily be applied to the wider range of meson states that we have considered. 

Recalling the ' = selection rule' [6, 12, 13], which strongly suppresses, by roughly 
two order of magnitude, the CEP cross section for non- = states, the Jz = vanishing 
of the gg — > MM production amplitude for flavour non-singlet scalar mesons will lead 
to a strong suppression in the predicted CEP rate. Moreover, more generally within the 
hard exclusive formalism, the production cross section of two mesons at high will be 
strongly suppressed by the numerically small factor {fM/k±)'^, where /m ~ 100 McV is 
the meson decay constant. These two combined effects have important phenomenological 
consequences, for example in the case of the potential tt^tt^ background to the candidate 
77 CEP events observed by CDF [5,16,20]: in particular, the possible tt^tt^ contribution 
is predicted to be very small. 

A further consequence of our results is that the non-resonant background to XcO ~^ '^'^ 
CEP within this perturbative framework is predicted to be small. However, in the re- 
gion M^TTT ~ it is not completely clear that the perturbative CEP mechanism should 
dominate. We therefore considered in Section 6 the possible 'non-perturbative' double 
Pomeron exchange contribution, calculated within the framework of Regge theory, which 
may give a comparable contribution at these Mj^^^ values, although there are some impor- 
tant uncertainties in the theoretical predictions at these experimentally unexplored regions 
of phase space. In principle, effects coming from the interference between the resonant 
Xc and continuum contributions should also be considered. The perturbative tail to this 
non-perturbative mechanism is also considered in Section 7. Although in general we find 
that it represents a power correction to the usual perturbative CEP process, which will be 
negligible in the limit of large meson k_L, in the case, for example, of the strongly suppressed 
TTTT CEP cross section, we find this may also be relevant at lower values of vrvr invariant 
mass. We shall consider these issues in more detail in a forthcoming publication [87]. 

Finally, we have also calculated the gg — > MM amplitudes for the production of flavour 
singlet meson states (e.g. r)'r]'), for which a new class of 'ladder' diagrams contributes, see 
Section 4, with the result that the Jz = amplitudes do not vanish. From this fact we 
predict that the CEP cross section for the production of scalar flavour singlet states should 
be strongly enhanced. This also applies to the case of r]r] and rji]' CEP, with the precise 
level of enhancement depending on the level of the t] — ij' mixing. A further interesting 
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possibility we have noted is that a valence gg component of the flavour singlet states may 
further increase the cross section. In this case, there is a large level of uncertainty in the 
precise size of such a component, and so the CEP mechanism could in principle shed light 
on this issue. 
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A. MHV formulae 

The spinor products are defined in terms of the solutions of the Dirac equation by 

{kikj) = {kr\kp =u_{ki)u+{kj) = v^{ki)v^{kj) , (A.l) 

[ki k,] = {k+\kr) = u+{ki)u-ikj) = v-{ki)v+{kj) , (A.2) 

with 

{ki kj)\kj kj\ = \ {ki kj)\ = 2(^kikj) . (A. 3) 

For more details and properties of the spinor products, see for example [61,95]. For com- 
pleteness we also reproduce the Schouten identity 

{ki k2) {ks ki) = {ki ki) (fcg ^2) + {ki ks) (fcs k^) , (A.4) 

B. Power-suppression of symmetric CEP mechanism: derivation. 
B.l Skewed mechanism 

The outgoing quark momenta are given by (see Fig. 15 (a)) 

li = xks l2 = {^-y)ki h = yki l4 = {l-x)k3, (B.l) 

where ^3^4 are the meson momenta. Considering the limit of forward outgoing protons, the 
incoming gluon momenta are given by 

qi = xipi + Q± q2 = X2P2 - Q± ■ (B.2) 
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We have 

= {xm + Q±- xksf ~ kl , (B.3) 

where the mesons have transverse momentum k±. Similarly we have k'^ ~ k'j_. The two 
quark propagators will therefore together give a factor of ~ 1/^1 to the amplitude. Also 

/c^ = t/(l-x)Mi~Ml, (B.4) 

where Mx is the meson pair invariant mass. Finally in the amplitude we have to include 
the polarization vector of the incoming gluons, 

2 2 

ei^t^2u -P1mP2i/ j^Q±nQ±u ~ , (B.5) 

where we have made use of the gauge invariance of the gg ^ X vertex q'^Vf^v = ^2^1^ = 0- 
Putting this together, the skewed CEP amplitude is given by 



fdQlQll I 1 1 Z-dQi 
B.2 Symmetric mechanism 

For forward outgoing protons, the fusing gluon momenta are given by (see Fig. 15 (b)) 

qi = xipi + gi_L qi = xipi - qi± , 

q2 = X2P2 + q2± h = X2P2 - q2± , (B.7) 

where gi,2 are defined as in Fig. 15 (b). As before we have k''^ ~ /c^ ~ A;^ and we pick up a 
factor of l/k'j_ from the quark propagators. qi± is unconstrained and has to be integrated 
over, while q2± is given by 

q2± = [x-{l- y)]k^ - q1 , (B.8) 

apart from near the constrained region of integration in the meson wavefunctions when 
X = {l — y)[i- + 0{qi±/k±)]. We therefore have ~ k'j_, i.e. for general quark 4-momenta 
there has to be a large momentum transfer through both gluons q2 and q2 (or of course qi 
and qi). From the gluon polarization vectors we pick up factors of qf 2±/-^x^ therefore 
the total amplitude is given by 



We therefore have 

^sym {Q±) 
^skew ^_|_ 

where 



(B.IO) 



iQl) - , (B.n) 
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with the usual Sudakov factor weights etc. implicit in the integrand. In general we have 
{Q\) = 3 — 4GeV^, depending on the central object mass Mx and cms energy -y/s. 

Finally, in the case of flavour-singlet meson production, for the additional diagram 
where the qq pairs forming the mesons are connected by a quark line, instead of (B.8) we 
have 

g2± = fc± - g± , (B.12) 

that is we have ~ over the full region of meson x, y integration, leading to a slightly 
larger suppression. 
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